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ABSTRACT 

We present optical/ncar-IR integral field unit (IFU) observations of a gas pillar in 
the Galactic Hn region NGC 6357 containing the young open star cluster Pismis 24. 
These observations have allowed us to examined in detail the gas conditions of the 
strong wind-clump interactions taking place on its surface. 

By accurately decomposing the Ha line profile, we identify the presence of a narrow 
(~20 kms -1 ) and broad (50-150 kms -1 ) component, both of which we can associate 
with the pillar and its surroundings. Furthermore, the broadest broad component 
widths are found in a region that follows the shape of the eastern pillar edge. These 
connections have allowed us to firmly associate the broad component with emission 
from ionized gas within turbulent mixing layers on the pillar's surface set up by the 
shear flows of the winds from the O stars in the cluster. We discuss the implications of 
our findings in terms of the broad emission line component that is increasingly found 
in extragalactic starburst environments. Although the broad line widths found here 
are narrower, we conclude that the mechanisms producing both must be the same. 
The difference in line widths may result from the lower total mechanical wind energy 
produced by the O stars in Pismis 24 compared to that from a typical young massive 
star cluster found in a starburst galaxy. 

The pillar's edge is also clearly defined by dense (<5000 cm~ 3 ), hot (>20 000 K), 
and excited (via the [Nil] /Ha and [Sll]/Ha ratios) gas conditions, implying the pres- 
ence of a D-type ionization front propagating into the pillar surface. Although there 
must be both photoevaporation outflows produced by the ionization front, and mass- 
loss through mechanical ablation, we see no evidence for any significant bulk gas 
motions on or around the pillar. We postulate that the evaporated/ablated gas must 
be rapidly heated before being entrained. 

Key words: ISM: H n regions - ISM: individual: NGC 6357 - ISM: kinematics and 
dynamics. 



1 INTRODUCTION 

1.1 Broad emission features 

Direct evidence of interactions between stellar winds 
and their surroundings is relatively sparse. Observa- 
tions of H2 and CO absorption lines on sight-lines 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile un- 
der programme 081.C-0808. 
f E-mail: msw@star.ucl.ac.uk 



towards massive stars in our Gal axy show compara- 
tively broad underlying components (jFalgarone fc Phillips! 
Il990l : iHartauist. Dyson fc Williams! 1 19921 ). These were in- 
terpreted as arising in turbulent surface layers on cold 
dense molecular clouds surrounding the stars result- 
ing from the impact of the stellar winds. Broad un- 
derlying components have also been observed in the 
optical emission lines from more en ergetic environ- 
ments in nearby giant Hn regions (e.g. Roy et all 1 19921 : 



IChu fc Kennicuttl Il994l; I Yang et al 
galaxies fe.g. iHomeier fc Gallagher 



1996) and starburst 



19991 : iMarlowe et all 
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1995; iMendez fc Estebanl 1 19971 ; | Sidoli. Smith fc Crowtherl 
20061 ) 



2006; Vanzi ct al. 



However, due to mismatches in 
spectral and spatial resolution and in the specific environ- 
ments observed, the nature of the energy source for these 
broad optical lines has been contested. Recent detailed in- 
tegral field unit (IFU) studies of the ionized interstellar 
medium (ISM) in the nearby starburst galaxies NGC 1569 
and M82, however, h ave shed a considerable am o unt of 
light o n this problem (|Westmoquette et al.1 l2007al lbl. |200S| . 
l2009al lbh. 

By mapping out the properties of the individual emis- 
sion line components in NGC 1569 (including a broad 
underlying component with FWHM ~ 100-350 kms -1 ), 
IWestmoquette et al.l (|2007al ) identified a number of corre- 
lations that allowed them to determine the likely origin of 
the broad component. They concluded that, like in the CO 
line case, the evaporation and/or ablation of material from 
cool interstellar gas clouds caused by the impa ct of the high- 
energy photons and fast-flowing cluster winds (|Pittard et al.l 



120051 ) produces turb ulent mixing layers (TML s ) on the sur- 
faces of the clouds (fBegelman fc Fabianlll990l ; ISlavin et al.l 
Il993l ; IBinette et all Il999l. 120091') fr o m whi ch the emission 
arises. IWestmoquette et al.1 (|2007d . l2009al lbh showed that 
this explanation is also applicable to the broad emission lines 
detected in M82. They argued that since M82's high pres- 
sure ISM is highly fragmented with many small clouds well 
mixed in with the star clusters, there are many cloud sur- 
faces with which the copious ionizing photons and fast winds 
can interact and this can explain the prominent broad line 
emission. 

IBinette et all (|2009l ) modelled the Ha and [Om]A5007 
broad emission components observed in the extragalactic gi- 
ant H 11 region NGC 2363 with a TML model. By using an 
approxim ate treatment of turbu lence (the mixing length ap- 
proach of ICanto fc Ragal Il99ll ). they avoided dealing with 
the details of the turbulent energy cascade, while still be- 
ing a ble to predict the l ine kinematics. The Mappings Ic 
code (|Ferruit et al.lfl997l ) was used to compute the radiative 
transfer and thus the line emissivities. In this model the 
line width is created by the velocity gradient in the mixing 
layer. A limitation of this approach is that the turbulent ve- 
locity that should be generated in the TML is not included 
and neither are the non-equilibrium ionization (NEI) effects 
that are expected when cooler gas mixes rapidly with hotter 
gas. This model is therefore only a good approximation if 
the mixing occurs fairly slowly and the ionization is domi- 
nated by photoionization. The model of lSlavin et al.l (|l993l ) 
includes NEI effects, but does not predict line profiles since 
the dynamics in the layer are not modelled. 

The only way to connect the low-energy molecular line 
cases to the higher-energy optical emission line cases is to 
model the physical process using the correct initial condi- 
tions. We have begun to do this using a hydrodynamical 
code that includes NEI effects, but we must first confirm 
that the physical mechanism producing the broad compo- 
nent in both cases is wind-driven TMLs, and we must have 
reliable measurements of the gas conditions within these lay- 
ers. Thus we have undertaken a study of a nearby Galactic 
analogue to the type of environments found in optical stud- 
ies of extragalactic giant H 11 regions and starburst galaxies. 
Although the energies and ISM conditions found in starburst 
galaxies are not easily found within our Galaxy, the advan- 



tage of proximity allows discrete superimposed kinematical 
components to be clearly resolved and unambiguously dis- 
entangled from other possible line broadening mechanisms. 

1.2 NGC 6357 

The analogue we have chosen is a gas pillar that forms 
part of the youngest Hn region (G353.2+0.9) in NGC 6357 
(RCW 131, W 22, Sh-2 11), a large Hn region/molecular 
cloud complex in the Sagittarius spiral arm of the Milky 
Way. G353.2+0.9 is ionized by the open cluster Pismis 24. 
This cluster has a stellar population a few times that of 
the Orion star cluster and a central stellar density sim- 
ilar t o that of the Gal actic young massive cluster NGC 
3603 (|Wang et al.1 120071 ). Pismis 24 is still too young t o 
have had any supernova explosions f Wang et al.l 120071 ); 
this is in-keeping with its O-star population, which in- 
cludes an 03.5 If (Pismis 24-1) an d 03.5 III(f*) (Pis- 
mis 24-17) star dWalborn et al.l I2OO2T ), each of MOO M 
(|Mafz Apellaniz et alj|2007l ). 

The conclusion of a number of individual radio and IR 
studies of NGC 6357 is that most of the molecular material 
is located either behind the Hll region and/or north of the 
sharp ionization front clearly visible in Fig. [T] (|Massi et al.l 
Il997h . The expansion of the Hn region has caused gas to 
accumulate in a (large) interim zone into which the ion- 
ization front is proceeding; this has triggered the forma- 
tion of several more (less massive) stars within the cloud 
jBohigas et al.ll2004 ). The molecular gas observations sug- 
gest that NGC 6357 represents a late-stage blister-type H 11 
complex viewed face-on, meaning that we are se eing the top 
ioniz ed face of a much larger molecular cloud (|Massi et al.l 
Il997l ). and that the pillar we have focussed on must project 
up from this structure towards the observer. Since the most 
highly excited regions of the nebula are facing the cluster, it 
is thought that the primary ionizing sources must be the two 
03.5 stars in the Pismis 24 cluster. This is confirmed by ion- 
ization calculations that show that t he exciting stars mus t 
be of spectral type earlier than 05 (|Bohigas et al.ll2004l ). 
The location of these stars and their projected distance to 
the gas pillar and ionization front (3.5 pc) are indicated in 

Fig.m 

The mechanisms that lead to the formation of "ele- 
pha nt trunks" in H 11 regions is still a matter of debate 
(e.g. iMizuta et al. I l2006l ; iGritschneder et al.ll2009h . However, 
the fact that these pillar-shaped nebulae usually point to- 
wards and exhibit velocity gradients in the direction of a 
nearby OB star, and show signs of compression and/or star 
formation at their tips, suggests a strong int eraction with 
the winds and radi a tion from these sources |Pound 1998 



Sugitani et"al"1l2002l; iGahm et alj|2006l; IMizuta et al 



2006 



Gritschneder et al.ll2009l ). Our targeted pillar in NGC 6357 



both points directly towards (in projection) the aforemen- 
tioned O stars in Pismis 24 (Fig. [TJ , and has been detected 
in both CO and HCO + impl ying that it contain s very high 
density gas (>10 5 -10 6 cm~ 3 ; lMassi et al.|[l997h . The pillar 
is also known to house two compact radio sources thought 
to be young stellar objects, including one at its apex (IRS 
4; IPersi et al.lll986l ; iFelli et al.lll99ol ; IWang et al.ll2007t ). 

In this paper we present optical/near-IR 
VLT/FL AMES- ARGUS IFU observations of the tip of 
this gas pillar (Section [2j . In these observations we have 
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Figure 1. Composite image of NGC 6357: the upper half comprises an ACS+WFPC2 F658N (Ha) composite. Since this image does 
not cover the stars of the cluster, we have extended the figure southwards by including part of the HST / ACS WF C F850W (~J-band) 
image. The ARGUS IFU position is outlined with a solid rectangle (11.5 X 7.3 arcsec), the two main ionizing stars JWalb orn et aLll2002| ; 
IBohigas et~al]|2004l) and the ionization front arc labelled, and the approximate projected distance from the stars to the tip of the pillar 
is indicated. Inset left, top: zoom-in showing the outline of the IFU position on the pillar. The image intensity scale has been stretched 
to enhance the protrusion indicated by the arrow and discussed in the text. Inset left, bottom: integrated Ha flux (C1+C2+C3) map as 
seen through the IFU. Ha line profiles and Gaussian fits from the two spaxels indicated are shown in the top right inset, where the flux 
units are in an arbitrary but relative scale. 



identified multiple broad and narrow emission line com- 
ponents. This has allowed us to examine in detail the 
kinematics (radial velocity and FWHM) of each Ha line 
component, and the spatially resolved nebular properties 
(densities, temperatures, excitations) of this region (Sec- 
tion [3]), and use these insights to discuss the wind and 
ionization interactions taking place on the pillar's surface 
(Section [4|. A summary of these findings and discussion is 
contained in Section [S] We adopt a distance to the NGC 
6357 complex of 2.56 kpc (|Massev et al.|[200ll ), meaning 10 
arcsec ~ 0.62 pc. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of NGC 6357 were obtained with the FLAMES 
instrument on the VLT with the ARGUS IFU coupled 
to the GIRAFFE spectrograph (PI: Westmoquette, Prop 
ID: 081.C-0808). We used the ARGUS array in its 0.52 
arcsec/spaxel spatial sampling mode, giving a field-of-view 
(FoV) of 11.5 x 7.3 arcsec sampled by 333 fibres, and centred 



it on 17 h 24 m 45=8, -34° ll'04'.'l (J2000) in order to cover 
the tip of the pillar. Image quality during the observations 
fell in the range 0.7-1.0". We observed this position with 
three grating settings (L614, L682 and L881; giving wave- 
length coverages of 5750-7200 A and 8200-9400 A) with ex- 
posure times of 2100-2700 sees. For each grating setting, two 
individual exposures were taken with a dither offset of 0'.'52 
in RA (1 fibre) in order to facilitate the masking of the dead 
fibres. Fig.Q]shows the position of the IFU on an HST ACS 
and WFPC2 F850W (~J-band) and F658N (Ha) composite 
image of NGC 6357, together with the location and distance 
of the pri mary ionizing source s of this nebula (Pismis 24-1 
and 24-17: IBohigas et al.|[2(5oi l. 

Basic data reduction was achieved with the girBLDRS 
(Giraffe BaseLine Data Reduction Software). This included 
bias subtraction, flat fielding, identification and spectrum 
extraction, and wavelength calibration. Post-processing was 
carried out in iraiQ, including cosmic-ray rejection (using 

1 The Image Reduction and Analysis Facility (iraf) is dis- 
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LACOSMlC; lvan Dokkuml200ll ). flux calibration and the com- 
bining of the individual exposures (dead fibres were masked 
and removed when the dithered exposures were combined). 

The FWHM spectral resolution, measured from averag- 
ing fits to isolated arc lines over all spaxels, was found to be 
0.55 A (L614 grating), 0.55 A (L682 grating), and 0.89 A 
(L881 grating). 



2.1 Emission line fitting 

The high S/N and spectral resolution of the data have al- 
lowed us to quantify the emission line profile shapes to a 
high degree of accuracy. In general, we find the emission 
lines to be composed of two bright, narrow components 
(FWHM ~ 10-30 kms" 1 ; hereafter CI and C3) overlaid 
on a fainter, broad component (FWHM ~ 40-150 kms ; 
hereafter C2). Follow i ng the methodology first employed by 
IWestmoquette et al.l (|2007af ) . we fitted multiple Gaussian 
profile models to each emission line using an iDL-based \ 2 fit- 
ting package called PAN (Peak ANalysis; Dimeo 2005). Each 
line in each of the 333 spaxels was fitted using a single-, 
double- and triple-Gaussian initial guess. Line fluxes were 
constrained to be positive and widths to be greater than the 
instrumental contribution to guard against spurious results. 
To fit the [Sn] doublet, for each component we constrained 
the wavelength difference between the two Gaussian mod- 
els to be equal to the laboratory difference, and FWHMs 
to equal one another. Multi-component fits were run sev- 
eral times with different initial guess configurations (widths 
and wavelengths) in order to account for the varied profile 
shapes, and the result with the lowest x 2 wa s retained. 

Determining, firstly, how many Gaussian components 
best fit an observed profile, and secondly, which fits are 
most appropriate, is not trivial. To decide how many Gaus- 
sian components best fit an observed profile (one, two or 
three), we employed the statistical F-test. This ascertains 
if the x 2 improvement obtained when increasing the num- 
ber of fit components is sta tistically significant or not (see 
IWestmoquette et al.l l2007al . for a fuller description of this 
procedure) . 

Experience has taught us that a number of additional, 
physically motivated tests are needed to firstly filter out 
well-fit but physically improbable results, and secondly to 
assess which Gaussian profile belongs to which physical line 
component. For a fit to be accepted, we set the criteria that 
the measured FWHM had to be greater than the associated 
error on the FWHM result (a common symptom of a bad fit) . 
Of a triple-Gaussian fit, we specified that the broadest com- 
ponent should be assigned to component 2 (C2), and after 
that, the redder to be component 1 (CI) and the bluer to be 
component 3 (C3). The reasons for this choice will become 
apparent in the following sections. With double-Gaussian 
fits, we assigned C2 to be the broader of the two components 
and CI to the other, regardless of their velocity difference. 
This consistent approach helped limit the confusion that 



tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 




2 4 6 8 10 12 



-25 -15 -5 5 15 25 

Radiol velocity {km s" 1 } 

Figure 2. Ha radial velocities in units of kms -1 (corrected 
to i>lsr) for the three line components. Contours represent the 
summed (C1+C2+C3) Ha flux as mapped in the bottom-left in- 
set to Fig. Q] 

might arise during analysis of the results where discontin- 
uous spatial regions might arise from incorrect component 
assignments. The top-right inset to Fig. [T] shows an exam- 
ple of a triple- and double-Gaussian Ha line profile extracted 
from the data cube, and illustrates how the components were 
assigned. 

3 EMISSION LINE MAPS 
3.1 Ha kinematics 

Maps of the Ha radial velocity (corrected to «LSR,kinematic) 
and FWHM (corrected for instrumental broadening) for the 
three identified line components are shown in Figs.[2]and [3] 
The radial velocities are shown relative to vlsr to a i d com - 
parison with the CO measurements of iMassi et al.l (|l997h . 
On each map, contours of the summed (Cl+C2+C3) Ha 
flux distribution are shown to locate the tip of the pillar 
and help guide the eye. 

These maps reveal a number of important characteris- 
tics about the ionized gas dynamics. South of the pillar the 
Ha line is split into a bright, redshifted, narrow component 
(CI), a fainter, narrow component (C3) with velocities rang- 
ing from +5 to —20 kms" 1 , and a broad, mostly blueshifted 
component (C2). The top panel of the right-hand inset to 
Fig- \H shows an example line profile from this region (Spec- 
trum 82). 

Ha C3 exhibits a velocity gradient from the south-west 
to the north and east of the field, ranging from —20 to 
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Figure 4. Average Ha line profiles and corresponding Gaussian fits from (a) the diagonal N-E to S-W strip containing the broadest 
broad line emission (total 55 spaxels) - C2 (broad component) FWHM = 87 kms ; (b) the remaining area (total 278 spaxels) - C2 
FWHM = 68 kms" 1 . 
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Figure 3. Ha FWHM in units of kms -1 (corrected for instru- 
mental broadening) for the three line components, CI, C2 and C3. 
Contours represent the summed Ha flux shown in the bottom-left 
inset to Fig. [T] 

+5 kms - . The CO velocity of gas surrounding the pillar 
(~0 kms" 1 ; iMassi et al.lll997M is in agreement with the red- 
dest of these velocities, suggesting that there is some corre- 
spondence between the molecular and ionized gas, although 
the bluest C3- and the majority of the C2-emitting ionized 
gas must certainly be flowing towards the observer at veloci- 
ties up to ~20 kms -1 relative to the CO gas. No evidence of 
the reds hifted CI compone nt can be seen in the CO velocity 
maps of lMassi et all (|l997ft . We caution, however, that the 
spatial resolution of the CO observations is ~40", meaning 
that a detailed comparison is meaningless. 

Moving from here onto the pillar surface, the bright 
redshifted component disappears and the remaining narrow 



component now becomes assigned to CI (hence the apparent 
discontinuity in the CI map). The bottom panel of the right- 
hand inset to Fig. [1] shows an example line profile from this 
region (Spectrum 189). Here the velocity of CI is consistent 
with the velocity of C3 in spaxels adjacent to this region, 
leading us to conclude that the component that was assigned 
to C3 in the spaxels surrounding the pillar is in fact the same 
component as that which becomes re-assigned to CI on the 
pillar (as the redshifted component is now obscured). The 
velocity of both the narrow (CI) and broad (C2) Ha com- 
ponents on the pillar are in very good agreement with that 
of the correspond ing CO velocities (wlsr ~ —5.5 kms -1 ; 
IMassi et al.lll997l ). indicating that the ionized and molec- 
ular gas in the pillar must be physically contiguous. The 
similarity of the broad component velocities with the bluer 
narrow component (C3 surrounding the pillar and CI on the 
pillar) suggests that C2 is indeed emitted by turbulent gas 
associated with the pillar. We discuss this further in Sec- 
tion n 

From Fig. [3] the width of the narrow components (CI 
and C3) appears very constant over the whole field-of-view 
(FoV). The mean FWHM of the two narrow components 
is 20 ± 7 kms -1 (where the quoted uncertainty represents 
the standard deviation on the sample and quantifies the de- 
gree of consistency). The broad component (C2; FWHM > 
40 kms -1 ) is identified in the majority of spaxels across 
the FoV, and rises to widths of ~150 kms -1 in placed]. 
The broadest Ha C2 lines (FWHM ~ 90-150 kms -1 ) are 
found along a diagonal strip extending from the north-east 
of the FoV, tracing the eastern edge of the pillar structure, 
along to the tip, then extending off towards the south-west. 
This coincidence lends further support to the association of 
the C2-emitting gas with the pillar surface. Close inspection 
of the HST Ha image reveals a very low surface brightness 
protrusion to the pillar tip that is coincident with the south- 
west extension to the diagonal strip of broad line emission. 
This is highlighted in the left-hand inset to Fig [1] and dis- 
cussed further in Section^ Fig.Ushows the average Ha line 
profile and the associated multi-component Gaussian fits for 
the 55 spaxels in the diagonal N-E to S-W strip containing 
the broadest broad line emission (panel a) and the remain- 



2 The broad component can also identified in the [Sm]A9069 line 
but at a lower signal-to-noise (S/N). 
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ing 278 spaxels (panel b). Although on average the broad 
component in this diagonal strip is fainter than in the rest 
of the field, it is clearly significantly broader (87 kms -1 vs. 
68 kms -1 ). 

3.2 Nebular properties 

The integrated (C1+C2+C3) Ha line flux map is shown in 
the bottom-left inset to Fig. [T] The vast majority of the Ha 
flux across the IFU field originates from CI and C3 (and 
then mostly CI); C2 makes a very small contribution to the 
total line flux [< 0.3 x (C1+C3)]. On the pillar, however, 
the narrow line (CI) is much fainter, and C2 becomes more 
dominant, but in some places C2 and CI make an almost 
equal contribution to the total line flux. 

Maps of the [N n]A6583/Ha and [S n]A6717+A6731)/Ha 
line ratios are shown in Figs. [5^. and [5p. Although nei- 
ther indicators show any evidence for shock-like ratios 
[log([Nll]/Ha)>-0.2 or log([S n]/Ha)>-0.4], both ratios 
are significantly enhanced on the pillar compared to the sur- 
rounding gas. In [Nll]/Ha (Fig. [5^), the peaks are located 
on the south-western and eastern pillar edges, whereas in 
[Sll]/Ha (Fig. [5p) the peak is located nearer the pillar's 
central axis, interior to the [Nil] /Ha peaks. This difference 
is made clear by comparing the [Sll]/Ha map to the over- 
plotted contours representing the [Nil] /Ha of Fig. [5^. Since 
the [S n] line flux may be enhanced by emission from par- 
tially ionized regions, we might expect [S n] to originate from 
deeper within the cloud compared to [N n] . The line ratio 
maps are certainly consistent with this idea. We might then 
also expect [S n] to exhibit a narrow line width, since the 
interior regions of the clump are presumably less turbulent. 
To investigate this we summed 20 spaxels covering the peak 
in [Nil] /Ha and 20 spaxels covering the [Sll]/Ha peak and 
measured the widths of the [N n] and [S n] lines in the two re- 
sulting spectra. To within the uncertainties, we find the line 
widths are equal and do not exhibit any significant change 
between the two locations. Observations with a S/N ratio 
sufficient to accurately decompose the [N n] and [S n] lines 
would be needed to investigate this in more detail. 

A map of the electron densities, calculated from the 
flux ratio of the [S ll]AA6717,6731 doublet (assuming T c = 
10 4 K), is shown in Fig. [SJ:. Here we sum the fluxes from 
all components in the [Sn] lines to increase the S/N ratio. 
The density of the gas projected on the body of the pil- 
lar is ~1000-1500 cm -3 , only marginally higher than that 
of the gas surrounding the pillar (500-1000 cm -3 ). These 
values are not unsurprisingly high for a large, gas-rich H II 
region. However, the measured densities at the pillar edges 
are significantly enhanced, peaking at ~5000 cm -3 . 

Electron temperatures were derived from the dered- 
dened flux ratio of [S m](A9069+A9532)/A6312. To dered- 
den this ratio, we first calculated the extinction in each 
spaxel using an average of the Pll/Ha, PlO/Ha, P9/Ha 
decrements, u sing the total summe d flux (C1+C2+C3) from 
each line, the ICardelli et al. ( 1981: ) extinction curve assum- 
ing Rv=3.5 ( Bohigas et al.l 2004 ). and the electron den- 
sities calculated above. We find extinctions in the range 
Ay — 3-5 mags, which i s consistent with previous studies 
(e.g. iBohigas et alj|2003 ). Since our observations only cover 
[S hi]A 6312 and [Sm]A9069, we a ssumed [S m] A9532/A9069 
= 2.48 (|Mendoza fc Zeippenll98l ) to recover the [S m]A9532 



fluxes. The resulting electron temperature map is shown 
in Fig. [5ji. Electron temperatures of ~12 000-14 000 K are 
found for the gas surrounding the pillar. The temperature 
then clearly increases towards the edges of pillar, peaking 
on the eastern edge at values of >20 000 K. The tempera- 
ture of the gas projected on the centre of the pillar (~9000- 
10 000 K) is lower than the surrounding gas. 

In summary, our excitation ([Sll]/Ha and [Nil/Ha), 
electron density ( [S n] doublet), and electron tempera- 
ture ([Siii][A9069+A9532]/A6312) diagnostics reveal that 
the gas on the pillar edge is clearly more highly excited, 
denser (~5000 cm -3 compared to ~1000 cm~ 3 ), and hotter 
(20 000 K vs. 12 000-14 000 K) than its surroundings. This 
suggests that the gas conditions change dramatically within 
the surface layers of the pillar. 



4 DISCUSSION 

The dense, hot, and highly excited conditions of the gas on 
the pillar's edge are consistent with what might be expected 
after a D-type ionization fron t has propagated into the pillar 
|Osterbrock fc Ferland|[2006T l. A D-type ionization front can 
occur when a supersonic ionization front runs into a dense 
neutral gas cloud. The increase in density causes the ion- 
ization front to slow down enough to allow the ionized gas 
behind the front time to move in response to the heating. 
The gas velocity quickly becomes supersonic, and the ram 
pressure of this material streaming away from the front to- 
wards the star (down the density gradient) causes a shock 
wave to break off from the ionization front and begin to com- 
press the gas ahead of the front. The shock front gradually 
weakens and the ionization front continues into the clump 
as a strong D-type front with a large density jump. 

The density and temperature increases we observe on 
the pillar surface are consistent with a radiative shock, as 
would be expected in such a dense environment. In this case 
the degree to which these properties increase is governed by 
the magnetic field strength that acts to limit compression 
and dominate the post-shock pressure. The material stream- 
ing away from the cloud, either as a result of the aforemen- 
tioned ram pressure or from the ionization and accompa- 
nying heating, forms what is known as a photoevaporation 
flow. 

Besides photoevaporation and TMLs, hydrodynamic 
models of wind-clump interactions also predict that mate- 
rial should be stripped (ablated) from the surface layers by 
the impact of the ste llar wind and incorporated (entrained) 
into the wind flow (|Hartquist et al.l 1 19861 ; iMarcolini et al.l 
2005). However, disregarding the redshifted narrow compo- 
nent presumed to originate from background gas, we see no 
evidence for any significant bulk gas motions on or around 
the pillar (such as split line profiles or diverging velocity gra- 
dients). The radial velocity of the narrow and broad com- 
ponents are furthermore consistent with that of the CO gas 
indicating that the ionized and molecular gas are physically 
contiguous. This lack of bulk motions could be for a num- 
ber of reasons: no evaporation/ablation and entrainment of 
material is taking place; the bulk gas motions are parallel 
to the pillar surface and therefore in the plane of the sky 
and not visible in the radial velocity maps; or the gas is 
rapidly heated to T > 10 4 K before being entrained into 
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Figure 5. (a) log([N n] A6583/Ha) ratio; contours represent the corresponding summed Ha fluxes (bottom-left inset to Fig. [TJ. (b) 
log([S ll]A6717+A6731)/Ha ratio. The contours represent the [Nll]/Ha map and are shown to highlight spatial difference between the 
peak [Sll]/Ha ratios and the peak [Nll]/Ha ratios, (c) electron density map; contours represent the summed Ha fluxes, (d) electron 
temperature map derived from the dereddened [S III] (A9069+A9532)/A6312 flux ratio. Each map has been calculated using summed fluxes 
from all identified line components (i.e. C1+C2+C3). 



the flow making it invisible in the optical. That no outflow 
and entrainment of material is taking place is highly un- 
likely given the clear evidence of a strong ionization front 
and the mechanical energy of the winds from the O stars 
that must be impinging on the clump surface. It is also un- 
likely that the gas motions are all sufficiently parallel to the 
pillar surface not to be visible, so we have to conclude that 
the evaporated/ablated gas must be rapidly heated before 
being entrained (i.e. the gas heating time is fast compared 
to outflow time). Thermal conduction-driven heating could 
help here since this process is known to be able to rapidly 
heat material so th at it quickly becomes part of the ambient 
flow (|Pittardll2007h . 

Could the large turbulent velocities we observe on the 
outer, optically-emitting surfaces of the pillar drive turbu- 
lence deeper into the pil lar core? In their m olecular line ob- 
servations of the pillar, iMassi et al.l (|l997l ) report CO line 
widths of 2-3 kms -1 . These are indeed highly suprather- 
mal for a molecular gas at their measured temperature of 
~30 K (where the sound speed would be ~0.5 kms -1 ). Fur- 
thermore, the reported line widths of HCO + , a species that 
traces higher density gas, are much narrower, suggesting 
that the the degree of turbulence must fall off rapidly with 
depth into the pillar. While it is conceivable that the large 
turbulent velocities that we observe on the outer surfaces of 
the pillar drive the turbulence of these deeper, denser layers, 
the large difference in velocities (100 kms -1 vs. 3 kms -1 ) 
implies that the situation may not be as straightforward as 
that. We intend to examine this question in the future with 
detailed hydrodynamical models. 



Contrary to expectation the region of broadest Ha C2 
line widths does not wrap around the western edge of the 
pillar, but continues from the pillar tip in the south-west 
direction, following the faint protrusion to the pillar (high- 
lighted in Fig. [TJ. If this region of broadest C2 emission 
traces the TML set up by the shear stellar wind flow, then 
the strongest wind-clump interaction must be taking place 
along this strip. It is interesting that the protrusion does 
not stand out in any of our nebular diagnostic maps; we at- 
tribute this to the fact that emission from the protrusion is 
too faint to be detected in any line but Ha in our spectra. 
Deeper observations would be needed to confirm this. 

What might this protrusion be? Two possibilities are: 
(1) it is a part of the original molecular gas cloud from which 
Pismis 24 formed that remains to be eroded or (2) that it has 
formed as a result of the hydrodynamics of the wind-clump 
interaction. If there were some outflow from the pillar as 
a result of ionization/heating, then in general it would flow 
away from the high pressure zone where the heating is occur- 
ring, along the sides of the clump and be swept downstream 
by the wind. A highly unstable situation would occur at the 
tip, where the pressure would push the outflow directly into 
the wind, and, for a time before the ram pressure of the wind 
overwhelms that of the outflow, could conceivably produce a 
protrusion such as we see. The geometry might be such that 
this protrusion is shielding the western edge of the main pil- 
lar body from the oncoming winds therefore preventing the 
development of a TML on its surface. 

An alternative idea is that the protrusion represen ts one 
half of a jet (or a monopolar jet; iHennev et alj |2002l) from 
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the IR and X-ray detected YSO known to i nhabit the apex of 
the pillar, IRS 4 ( likely a B0-B2 protostar; [persi et al.|[l98l 
I Wang et all 120071 ): in this case the turbulent line widths 
would arise from interactions between the jet and the sur- 
rounding gas. Although the length of the protrusion (~2. 6 " 
« 3300 AU) is plausible for a microjet (|Ballv et aljbood ), 
the fact that the structure appears curved and that we see no 
evidence for it in the radial velocity maps makes this argu- 
ment very unlikely. We therefore favour the former scenario 
where the structure is formed from a remnant overdensity 
in the ambient medium or an instability in the outflowing 
material. 

The existence of a broad component, not just on the 
pillar surface, but across the entire IFU field-of-view indi- 
cates that turbulent 10 4 K gas pervades the entire region. 
Although detailed hydrodynamic models would be needed 
to ascertain whether this gas exists as a result of the evapo- 
ration/ablation flows from the pillar surface, the significant 
difference in width between the broad component associ- 
ated with the TML (>100 kms -1 ) and the surrounding gas 
(~50 kms -1 ), together with the excitation and density re- 
sults, suggests that these two gas bodies are not related. Can 
these results help us interpret what w e see in more distant 
and e ne rgetic starburst systems (e.g. IWestmoquette et al.l 
l2007al fbl. [2009aH bh? The C2/C1 flux ratio is, in general, sim- 
ilar to what we have found in starburst galaxies, suggesting 
that the TMLs in both cases are of similar importance. The 
C2 widths, however, are narrower (50-150 kms" 1 vs. 150- 
350 kms -1 ). The de gree of turbulence in the TML might be 
expected to depend on the ratio of mass flux in the wind to 
that in the ionization initiated outflow from the clump. We 
know the Pismis 2 4 cluster contains only a few O stars [Pis- 
mis 24-1 - which iMafz Apellaniz et al.l (|2007h resolve into 
three objects with a combined mass of ~200 M© - and Pis- 
mis 24-17] an d is too young to h ave had any supernovae (i.e. 
age <4 Mvr; Iwkng et ai1l2007l ). Assuming a mass loss rate 
of 2 x 10 -5 M© yr -1 and a terminal velo city of 2500 kms -1 
for an 03-type star l|Smith et alj|2002h . this gives a wind 
power of 4 x 10 37 erg s -1 per star or a total of 1.6 x 10 38 for 
4 stars. Starburst99 jLeitherer et al.lfl999l ) predicts a wind 
power of 10 40 erg s -1 for a 10 6 Mq super star cluster, such 
as those typically found in starbursts. The total mechanical 
wind energy emitted from the Pismis 24 cluster is thus at 
least two orders of magnitude less than this, meaning that 
a lower level of turbulence in the surface layers of the pillar 
compared to those found in starburst environments is not 
unexpected. 



5 CONCLUSIONS 

In this study we have examined in detail the gas conditions 
in the surface interaction layers of a gas pillar in NGC 6357 
ionized by the star cluster Pismis 24 using high-resolution 
optical/near-IR IFU observations. Our observations show 
that the gas on the pillar's edge is denser, hotter, and more 
highly excited than its surroundings, providing clear evi- 
dence that the incoming radiation and winds from the clus- 
ter stars are strongly affecting the state of the ionized gas 
in these surface layers. 

That we have been able to accurately decompose the 
Ha line into multiple components has allowed us to exam- 



ine the gas dynamics of this region in great detail. One of 
the primary aims of this study was to look for the presence 
of a broad emission component, and to investigate whether 
this component originates in turbulent mixing layers (TMLs 
ISlavin et ahll 19931 ; iBinette et aHl2009h on the surface of the 
pillar set up by the shear flows of the winds from the young 
massive stars of the cluster Pismis 24, s i nce th is is the 
scenario suggesed by [Westmoquette et all (|2007al , see also 
IWestmoquette et al.ll2007bl la. l2009al lbT) to explain the broad 
emission components found in the starburst galaxies NGC 
1569 and M82. On the pillar we find clear evidence of both 
a narrow (~20 kms -1 ) and broad (50-150 kms -1 ) com- 
ponent; away from the pillar we also see a much brighter, 
redshifted narrow component. The similarity of the broad 
C2 component velocities to the velocity of the narrow line 
component emitted at the location of the pillar leads us to 
conclude that the turbulent, C2-emitting gas must be as- 
sociated with the pillar structure, whereas the redshifted 
narrow component originates in the background gas. That 
the velocities of the components associated with the pillar 
are consistent with those of the CO gas, and the fact that 
the broadest C2 line widths are found to follow the shape of 
the eastern pillar edge, lends further support to the connec- 
tion between C2 and turbulent ionized gas within the pillar 
surface layers. We can therefore assert that the mechanical 
energy of the stellar winds impacting on the pillar surface 
is driving a turbulent mixing layer, and it is from this layer 
that the broad component is being emitted. 

The pillar's edge is also clearly defined by dense 
(<5000 cm -3 ), hot (>20000 K), and excited (via the 
[Nii]/Hq and [Sn]/Ha ratios) gas conditions, implying the 
presence of a D-type ionization front propagating into the 
pillar surface. Although there must be both photoevapora- 
tion outflows produced by the ionization front, and mass- loss 
through mechanical ablation (stripping) , we see no evidence 
for any significant bulk gas motions on or around the pil- 
lar. We postulate that the evaporated/ablated gas must be 
rapidly heated before being entrained (i.e. the gas heating 
time is fast compared to outflow time). 

Although the location of the broadest Ha C2 emission 
follows the shape of the eastern pillar edge, it continues 
from the southern tip along a diagonal strip extending in 
the south-west direction. This extension is coincident with a 
faint protrusion to the tip of the pillar seen in the HST / ACS 
F658N image. We conclude that it is either a part of the 
original molecular cloud that remains to be eroded or that 
it has been formed through a hydrodynamical instability in 
the outflowing material, and the geometry is such that it is 
shielding the western edge of the main pillar body from the 
oncoming winds therefore preventing the development of a 
TML on its surface. 

Finally we have discussed the implications of our find- 
ings to what we have previously observed in starburst galax- 
ies. The C2/C1 flux ratio is, in general, similar to what we 
have found in these energetic extragalactic environments, 
suggesting that the TMLs in both cases are of similar im- 
portance. That the C2 widths are narrower (50-150 kms 
vs. 150-350 kms -1 ), however, could result from the fact that 
the Pismis 24 cluster contains only a few O stars and is too 
young to have had any supernovae, so the total mass flux 
and therefore mechanical wind energy is at least two orders 
of magnitude lower than that from a typical young massive 
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star cluster such as those found in starburst galaxies. This 
effect results in a lower level of turbulence in the surface 
layers of the pillar. 

To further investigate the relationships between the dif- 
ferent gas components within the pillar, its surface and the 
surrounding gas more fully, it would be important to de- 
rive excitation, density and temperature diagnostics for each 
of the individual line components (CI, C2 and C3). This 
we have not been able to do since our observations of the 
faintest lines required for these diagnostics were not of a suf- 
ficient S/N. Thus, we intend to obtain deeper observations 
of this object in the future to build on the work presented 
here. 
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